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Abstract 

In this paper we present an algorithm to compute aU unirational fields 
of transcendence degree one containing a given finite set of multivariate 
rational functions. In particular, we provide an algorithm to decompose 
a multivariate rational function / of the form / = g{h), where <? is a 
univariate rational function and h a multivariate one. 

1 Introduction 

Let K be an arbitrary field and K(x) = K(a;i, . . . ,a;„) be the rational function 
field in the variables x = {xi, . . . , a;„). A unirational field over IK is an inter- 
mediate field F between K and K(x). We know that any unirational field is 
finitely generated over K (see [B]). In the following whenever we talk about 
"computing an intermediate field" we mean that such finite set of generators is 
to be calculated. The problem of finding unirational fields is a classical one. In 
this paper we are looking for unirational fields F over K of transcendence degree 
one over K, tr.deg(F/K) = 1. 

In the univariate case, the problem can be stated as follows: given univariate 
rational functions /i , . . . , e , we wish to know if there exists a proper 
intermediate field F such that K(/i, . . . , /„) C F C K(y); and in the affirmative 
case, to compute it. By the classical Liiroth theorem (see [M]) the problem 
is divided in two parts: first to compute / such that IK(/i, . . . , fm) = 1K(/), 
and second to decompose the rational function /, i.e., to find g,h G 1K(?/) such 
that, F = K{h) with / = g{h). Constructive proofs of Liiroth's theorem can be 
found in [7], |12| and [1]. Algorithms for decomposition of univariate rational 
functions can be found in [17j and [1]. 

In the multivariate case, the problem is: given /i, . . . , /,„ in ]K(x) we wish 
to know if there exists a proper intermediate field F such that IK(/i, . . . , fm) C 
F C IK(x) with tr.deg(F/K) = 1; and in the affirmative case, to compute it. A 
central result is the following generalization of Liiroth's theorem: 

Theorem 1 (Extended Liiroth's Theorem). Let ¥ be a field such that K C F C 
K(x). //tr.dcg(F/]K) = 1, then there exists f e K(x) such that F = IK(/). 
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Such an / is called a Liiroth's generator of the field F. This theorem was first 
proved in [2] for characteristic zero and in [4] in general, see also [11] Theorem 
3. Using Grobner basis computation, the paper [S] provides an algorithm to 
compute a Liiroth's generator, if it exits. See also [S] for another algorithmic 
proof of this result. In this paper we present a new algorithm, which only 
requires to compute gcd's, to detect if a unirational field has transcendence 
degree 1 and, in the affirmative case, to compute a Liiroth's generator. We also 
present a constructive proof of the above theorem for polynomials (see [5]): if 
the unirational field contains a non-constant polynomial, then it is generated by 
a polynomial. 

By the Extended Liiroth's theorem, to find an intermediate field of tran- 
scendence degree one is equivalent to the following: first to find a Liiroth's 
generator /, i.e., K(/i, . . . ,/,„) = ]K(/), if it exists, and second to decompose 
the multivariate rational function /, i.e., to find g G IK(y) and h S IK(x) such 
that / = g{h) in a nontrivial way. The pair (g, h) is called a uni-multivariate 
decomposition of /. We present two algorithms to compute a nontrivial uni- 
multivariate decomposition of a multivariate rational function, if it exits. 

This paper is divided in four sections. In section 2 we state the proof of the 
Extended Liiroth's theorem and its polynomial version. In section 3 we present 
and analyze two algorithms to compute a uni-multivariate decomposition of a 
rational function, if it exists. In section 4 we discuss the performance of these 
algorithms. 

2 The Extended Liiroth Theorem 

In this section we present an algorithm to the following computational problem: 
given fi, . . . ,fjn £ IK(x), to compute a Liiroth generator / for F = K(/i, . . . , fm) 
if it exists, moreover we detect if F contains a non-constant polynomial, and in 
the affirmative case we find a generating polynomial. We start with the following 
definition: 

Definition 1. Let p G IK[a;i, . . . , Xn, yi, . . . , y„] = K[x, y] be a non-constant 
polynomial. We say that p is near— separated if there exist non-constant poly- 
nomials ri, si G K[x] ~ K.[xi, . . . , Xn] and r2, S2 G K[y] = K[yi, . . . , ?/„], such 
that neither ri, si are associated, norr2, S2 cire associated andp — riS2—r2Si. In 
the particular case whenp = r{xi, . . . , Xn)s{yi, . . . , y„)-s(a:i, . . . ,x„)r(?/i, . . . ,?;„), 
we say that p is a symmetric near— separated polynomial. We say that (r, s) 
is a symmetric near separated representation of p. 

In this paper, deg^.^ will denote the total degree with respect to the 
variables xi, . . . ,Xn and deg will denote the total degree with respect to all the 
variables. Also, given a rational function / we will also de note as fm fd the 
numerator and denominator of /, respectively. 

In the following theorem we give some basic properties of near-separated 
polynomials, for later use. 
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Theorem 2. Let p G K[x, y] be a near-separated polynomial and ri,si,r2,S2 

as in the above definition. Then 

(i) //gcd(ri, si) = 1 and gcd(r2, .S2) = 1, p has no factors in K[x] or K[y]. 

(ii) deg^.^ = niax{deg ri,deg si} anddegy^ ^^p = max{deg r2,deg 52}- 

(iii) // p is symmetric and (ai, . . . , «„) G IK" verifies 
p(xi, . . . , Xn, ai, . . . , a„) ^ 0, then there exists (r, s), a symmetric near- 
separated representation of p, such that 

r(Qi, . . . ,a„) = and s(q;i, . . . , a„) = 1. 

(iv) If p is symmetric, the coefficient ofx]°y-j^ in p is the near-separated poly- 
nomial 

where ai is the coefficient of x\ in r and hi is the coefficient of x\ in s. 

Proof, (i) Suppose v G K[xi, . . . , Xn] is a non-constant factor of p. Then there 

exists i such that dcg^.v > 1. Without loss of gcncrahty wc will suppose that 
i = 1. Let a be a root of v, considering p as a univariate polynomial in the 
variable xi, in a suitable extension of lK[x2, . . . , x„]. If a is a root of any of the 
polynomials ri or si, then it is also a root of the other. This is a contradiction, 
because gcd(ri, si) = 1. Therefore a is neither a root of n nor of si. Then, 

ri{a,X2,.-.,Xn) ^ r2{yi, ■ ■ ■ ^ ^ 
si{a,X2, . . . ,Xn) 32(2/1, y„) 

A contradiction again, since r2 , S2 are not associated in K. 

(ii) If deg ri 7^ deg si, the equality is trivial. Otherwise, if deg n = deg si > 

deg^^ 2;„Pi the terms with greatest degree with respect to xi,. . . ,Xn vanish. 
This is a contradiction, because r2, S2 are not associated. The proof is the same 
for r2, S2. 

(iii) Let (r, s) be a representation of p. 

— If r(ai, . . . , ttn) = 0, since . . . , a;„, qi, . . . , a„) 7^ 0, we have s(ai, . . . , «„) 7^ 
0. Then we have a new near-separated representation: 

rs{ai, . . . ,q;„). 



s{ai, . . . ,a„) 

— If s(ai, . . . ,«„) = 0, then the representation (— s,r) we are in the previous 

case. 

— If r(ai, . . . , a„), s{ai, . . . , an) 0, then we consider the representation 

rs{ai, . . . ,Q!„) - sr{ai, . . . ,a„), — r 

s[ai,. . . ,a„j 
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(iv) This is a simple routine confirmation. 



□ 



Now, we state an important theorem that relates uni-multivariate decom- 
positions to near-separated polynomials, that is proved in [10] : 

Theorem 3. Let A = IK(x) and M = K(y) be rational function fields over 
K. Let /, /i G A and f',h' £ M be non-constant rational functions. Then the 
following statements are equivalent: 

A) There exists a rational function g G K(t) satisfying f = g{h) and f — g{h'). 

B) h- h' divides f - f m A (g)K B. 

As a consequence, a rational function / G K(x) verifies / = g{h) for some 
g, h if and only if /i„(x) hd{y) - /id(x) /i„(y) divides /„(x) /^(y) - /d(x) /„(y). 

Given an admissible monomial ordering > in a polynomial ring and a nonzero 
polynomial G in that ring, we denote by Im G the leading monomial of G with 
respect to > and Ic G its leading coefficient. 

Algorithm 1. 

Input: /i, . . . ,/„, G K(x). 

Output: / G K(x) such that K{f ) = F = K(/i, . . . , /„), if it exists. Otherwise, 
returns null. 

A Let > be a graded lexicographical ordering for 
y = (yi, • ■ • Let i = m. 

B Let Fk = fkniy) - /fe(x)A.rf(y) for k ^ 1, . . . ,i. 

C Compute Hi — gcd({i^fc, k = I, . . . , i}) with Ic Hi — 1. 

D — If Hi = 1, RETURN NULL (¥ does not have transcendence degree 1 
over IKJ. 

- // there exists j G {1, . . . , i} such that Im Hi = Im Fj , then RETURN 

h- 

— Otherwise, let fi+i be a coefficient of Hi in¥\ K. Increase i and go to 

B. 

Correctness proof. If F has transcendence degree 1 over K, we can write 
F = K(/). By Theorem H 

/n(y) - /(x)/d(y) 

divides Hi for any i. Therefore, Hi is non-constant if a Liiroth generator exists. 

If there are i,j such that Im Hi ~ Im Fj, then Fj is a greatest common 
divisor of {Fk,k = !,...,«}. Therefore, Fj divides F^ for every k. Fix such 
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a fc. Let q — fkniy) I '5 — fkdiy) the normal form with respect to the 
monomial ordering >; then there exist p, g, r, s G F[y] such that 

where Im Fj does not divide any monomial of g or s. By theorem [5{i), g, s ^ 0. 
By the definition in step B, 

=f;(p-r/fe(x))-(g-s/fe(x)). 

Hence Fj divides q — s fk{xi, . . . ,Xn) and we conclude that q — s/fe(x) = 0, 
since otherwise we would get that Im Fj divides Im [q ~ s fk{xi, . . . , a;„)), which 

contradicts the choice of the polynomials g, s. Thus fk{xi, . . . , Xn) = - e F = 

K(/,). 

Now we suppose that Im Hi < Im Fk for all k. Again, fix a value for k. Then 
there exists a C G F[y] \F such that Fk = HiC. Let d, a be the lowest common 
multiples of the denominators of the coefficients of Hi and C, respectively. Then 
D — Hid, C = aC G K[x, y]. Since Hi is monic, the polynomial D is primitive. 
Then, 

D C 

fkniy) fkdi^) - /feriW fkdiy) ^ -7 — fkd 

a a 

and by theorem [21 

fkniy) fkdi^) - A„(x) fkdiy) = DC, 

C e K[x, y]. On one hand, D ^ K[y], thus D (and Hi) have a non-constant 
coefficient. On the other hand, C ^ K[y], then the non-constant coefficients of 
D in the ring K(x)[y] have smaller degree than that of /fe(x). The choice of d 
assures that the coefficients of H have smaller degree than fk- Summarizing, 
there exists a coefficient a G F of Hi that can be added to the list of generators 
and has smaller degree than them. If tr.deg (F/K) = 1, is non-constant 
for all i, and the generator has smaller degree than the others. Therefore, the 
algorithm ends in a finite number of steps. 

Finally, we note that complexity is dominated in the step C by computing 
gcd's of multivariate polynomials, so the algorithm is polynomial in the degree 
of the rational functions and in n (see |16)). 

From the fact that the Liiroth generator can be found with only some gcd 
computations, we obtain that if / is a Liiroth generator of K(/i, ...,/„) then 
it is also a Liiroth generator of ]K'(/i, . . . , /„) for any field extension K' of K, 
K C K'. 

Example 1. Let Q{fi, f2) C Q{x,y,z) where 

y'^x^ — 2y^x^z + y'^z'^ + — 2xz + z^ 



fi 

f2 = 



yx^ — yxz — yzx^ + z^y 
y'^x^ - 2y^x'^z + y'^z^ 



2xz + yx^ — yxz + z^ — yzx^ + z^y 
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Let 

Ft = ftn{s,t,u) - fi{x,y,z)f^^{s,t,u) , I = 1,2. 

Compute 

u A^T? T?\ . , 2. x'^y-zy , -x^y + zy 

H2 = gcd(Fi, F2) = -tu + s i H u H s. 

X — z X — z 

Since Im H2 < Im Fi with respect to the lexicographical ordering s > t > u, we 

X y — zy 

take a non-constant coefficient of H2 ■' /s = Now 

X — z 

2 ^"^y ~ ^y "X^y + ^y 

ii3 — —tu + s < H u H s 

and = since = _ff2- The algorithm returns /a, a Liiroth generator of 

Q(/l,/2). 

It is important to highlight that when the field F contains a non-constant 
polynomial you can compute a polynomial as a generator, and this generator 
neither depends on the ground field K. This result was proved in [5?, for zero 
characteristic. A general proof can be found in [llj . 

Algorithm 2. 

Input: /!,...,/,„£ lK(x). 

Output: / G ]K[x] such that K(/) = F = K(/i, . . . , if it exists. Otherwise, 
returns null. 

A Compute a Liiroth generator f o/K(/i, . . . , /„) using Algorithm{J[ 
B Let s be the degree of f . 

— If s > deg /4 and f'^ is not constant, return null. Otherwise, let 

f^l/f- 

— If s > deg and f'^ is not constant, return null. Otherwise, let 

f^f- 

— Let fn"^ T fct'^ ^6 homogeneous components of degree s of fn^fd, 

fn^' 

respectively. Let a = . If a or f^ ~ af'^ are not constant, return 

■Id 

null. 

Otherwise, let f — o /'. 

y~ a 

Correctness proof. Once a Liiroth's generator has been computed, take a 
generator / with degree m such that if / = ^ and 

fd 
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f - &^ + • • • + f'-"^ 

Jn — Jn T I Jn J 
Jd -Id ^ ^ Jd ' 

the sum in homogeneous polynomials, then either = or — ^ ^ K. 

If p e ]K(/) is a polynomial, then there exists g G ^{y) with degree r such that 
o-rV^ + ■ • • + Oo 



p = 5(/). If 5 = 



P 



Wif + ■ • • + 60 ' 

Or/; + • • • + ao/d 



&r /;: + ••• + foo/,^ 
Kifk^ + ■■■ + fi°^Y + ■■■ + b,if^/^ + ■■■ + ff^y ■ 

Since p is a polynomial, the degree of the previous denominator is smaller than 

the degree of the numerator. Therefore brfn + • • • + &o/d — 0- 

If ^ then 6^ = and p = + ' ' ' + ""-^^ Hence fd divides 

' /d(6.-i/r^ + --- + &o/r ) 

the numerator of p, and therefore divides /„. This proves that / is a polynomial. 

If, on the contrary, f^' 7^ 0, (jd I —j^ 1 = 0. Contradiction, since K. 

\fd / fd 



3 Two uni— multivariate decomposition algorithms 

We define the degree of a rational function / = fn/ fd G IK(x) as deg / = 
max {deg /„, deg fd} if gcd(/„, fd) = 1. The following definition was introduced 
in 15J for polynomials. 

Definition 2. Let f,h <E IK(x) and g G K(2/) such that f = g{h). Then we 
say that {g,h) is a uni-multivariate decomposition of f . It is non-trivial if 
1 < deg h < deg /. The rational function is uni-multivariate decomposable if 
there exits a non-trivial decomposition. 

If / is a polynomial having a nontrivial uni-multivariate decomposition, then 
by Algorithm [5] we get that there exits a uni-multivariate decomposition {g, h) 
with g and h polynomials. The paper [15] provides an algorithm to compute a 
nontrivial uni-multivariate decomposition of a polynomial / of degree m that 
only requires 0{nm{m -\- l)"log m) arithmetic operations in the ground field 
K. 

The known techniques for decomposition all divide the problem into two 
parts. Given /, in order to find a decomposition / — g{h), 

1. one first computes candidates h, 
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2. then computes g given h. 

Determining g from / and /i is a subfield membership problem. The paper 
[13] gives a sohition to this part. We also present another faster method, that 
only requires solving a linear system of equations. Usually, the harder step is 
to find candidates for h. One goal in decomposition is to have components of 
smaller degree than the composed polynomial. This will be the case here. 



3.1 Preliminary results 

First, we state some results that will be used in the algorithms presented later. 
On the properties to highlight out of uni-multivariate decomposition is the good 
behaviour of the degree with respect to this composition. 

Theorem 4. Let g e K(y) and h e IK(x), and f = g{h). Then dog / = 
deg g ■ dcg h. 

Proof. Let g = — and h = with gcd((7„, g^) = 1 and gcd{h„ , hd) = 1. Then 

gd hd 
there exist polynomials A, B £ K[y] such that 

9n{y) A{y) + gaiy) B{y) = 1. 

Homogenizing the polynomials g„, g^, A, B we obtain, respectively, the bivariate 
polynomials gn{yi,y2),gd.{yi,y2), A{yi,y2),B{yi,y2) verifying 

gn{yi,y2) A{yi,y2) y2+9d{yi,y2) B{yi,y2) yl = 
with either u = or v = and w =max{u, w}. Therefore, 

gn{hn, hd) A{hn, hd) hd +gd{hn, hd) B{hn, hd) hd = hd- 

If d is an irreducible factor of gcd(g„(/i„, hd),'gd{hm hd)), then d divides hd- On 
the other hand, d divides gn{hm hd) and gd{hn, hd)', this implies that d divides 
/i„. As a consequence, gcd(g„(/i„, /i^), /i-d)) = 1- So, 

f^ 9nf-M) |a| = |deg/.„-deg/.,| 
gd{hn, hd) 

is in reduced form. Without loss of generality, we can take deg gn = fn > 
Td =deg gd with 

9n{y) = y''" H h ao 

gd{y) = W^y''^ + ■■■ + bo. 

Then, deg / = max {deg 5„(/i„, hd), deg gdihn, hd) /i^" '^''l ^'^'^ 

gn{hn, hd) = H h aoh^i" 

gd{hn,hd) ^ hr^Ki^ H ^boh'2'. 
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If deg gn{hn, hd) = rn dcg h, we immediately obtain that deg / = deg gdeg h. 
If deg gn{hn, hd) < r„deg h, then s = deg /i„ = deg hd- Write 

/»„ = hi''' + h^r'^ + • • • + /ii"' 
hd^h^^ + h^-^ + --- + h^ 

(i) (?) 

where /i„ , hj are the homogeneous components of /i„ and /id with degree j, re- 
spectively. Since the degree decreases, deg g„(/il''\ ft, |,*^) = Oanddegg„ ( — ^ 1 ~ 

hi') 

0. Therefore, — ^ £ K. In this case, you can take h' G K(x) a rational func- 

^d 

tion with deg h = degh' , deg h'^ ^ deg h'^ and such that / = g'ih!) for some 
g' G with deg g = deg g'. Under these hypothesis, we proved before that 
deg / = deg g' deg h' — deg 3 deg h. □ 

Corollary 1. Let g — gnl gd with gn — a„j/" + • ■ • + ag, gd = by-y^ + • ■ ■ + &o o.i^'d- 
h = hn/hd verifying gcd{gn, gd) = gcd(/i„, hd) = 1. // / fn/fd = g{h) with 

f,, = {ayhl + --- + a,K)hT''^^-''^'^ 
fd = {byhl + --- + b,hl)hT''^-'''''^ 

then gcd(/„, /(;) = !• 

Proof. It is easy to prove that 

deg /„, deg fd < max{u, v} ■ max{deg ft„, deg hd}. 

If gcd(/„, fd) ^ 1, then dcg / < deg gdeg h, contradicting theoremHl □ 

Corollary 2. Given f and h, if there exists g such that f — g{h), then g is 
unique. Furthermore, it can be computed from f and h by solving a linear system 
of equations. 

Proof. If / = gi{h) = g2{h), then {gi - g2){h) = 0, and by theoremH deg (51 - 
52) — 0, thus gi — g2 is constant. It is then clear that it must be 0, that is, 
gi ~ g2. Again by theorem 31 the degree of g is determined by those of / and h. 
We can write g as a function with the corresponding degree and undetermined 
coefficients. Equating to zero the coefficients of the numerator of / — g{h), we 
obtain a linear homogeneous system of equations in the coefficients of g. If we 
compute a non-trivial solution to this system, we find g. □ 

Definition 3. Let f G K(x) be a rational function. Two uni-multivariate de- 
compositions {g, h) and {g' , h') of f are equivalent if there exists a composition 
unit I G ]K(?/), i.e., deg I = 1, such that h = l{h'). 

Corollary 3. Let f G IK(x) be a non-constant rational function. Then the 
equivalence classes of uni~multivariate decompositions of f correspond bijec- 
tively to intermediate fields ¥, K(/) C F C ]K(x), with transcendence degree 1 
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Proof. The bijection is 



{[i9,h)], .f = gm 
[i9,h)] 



{K(/) c F, (F/K) = 1}. 
F = K{h) 



Suppose we have a uni-multivariate decomposition (<?, /i) of /. Since / = g{h), 
¥ = K(/i) is an intermediate field of K(/) C IK(x) with transcendence degree 1 
over K. Also, if {g', h') is equivalent to (g, ft-), h = l{h') for some composition unit 
I e K(y). Consequently, h' = I'^ih) and ]K(/i) = ]K(/i'). Let (g, /i) and {g',h') 
be two uni-multivariate decompositions of / such that = K(/i'). Then 

there exists rational functions 1,1' G lK(y) such that /i = l{h') and ft' = 
By theorem m deg l{l') — 1 and deg I — deg V = 1. By the uniqueness (see 
Corollary [2]) of the left component, y = So, I G IK(y) is a composition 

unit and (g, ft), ((/', ft') are equivalent. By Theorem[TJ there exist ft G IK(x) and 
g e ]K(?;) such that F = K(ft) and / = .g(ft). □ 

3.2 First algorithm 

We now proceed with the first algorithm for computing candidates ft — hn/hd- 
This algorithm is based on Theorem [31 Since ft„(x) ftd(y) — ftd(x) ftn(y) divides 
/„(x) /d(y) — /d(x) /n(y), one can compute candidates for h from / merely 
looking at the near-separated divisors H = r(x) s(y) — r(y) s(x). Next, the 
problem is: given a multivariate polynomial H ~ (x, y) , how can one determine 
if it is a symmetric near-separated polynomial? This is a consequence of theorem 



Corollary 4. Given a polynomial p G IK[x, y], it is possible to find a near- 
separated representation (r, s) G K[x]^ of p, if it exists, by solving a linear 
system of equations with coefficients in K. Moreover, any other solution (r', s') 
of this linear system of equations gives an equivalent decomposition. 

Algorithm 3. 

Input: / G K(x). 

Output: A uni-multivariate decomposition {g,h) of f, if it exists. 
A Factor the symmetric polynomial 



Let D = {Hi, . . . , Hm} the set of factors of p (up to product by constants) . 
Let 1 = 1. 

B Check if Hi is a symmetric near-separated polynomial. If H = r(x) s(y) — 

r 

r(y)s(x), then h = — is a right-component for f ; compute the left compo- 
s 

nent g by solving a linear system (see Corollary\^ and RETURN (g, ft). 



m 



P = /«(x) /rf(y) - /<i(x) /„(y). 
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C Ifi<m, then increase i and go to B. Otherwise, RETURN NULL (f is 

uni-multivariate indecomposable ) . 

Example 2. Let 

f = Az^y'^ — 8z^y^ + Sz^yx + iz^y^ — 8zy^x 
+4a;2 - 2z'^y + 2zy'^ - 2x + 10. 

The factorization of the polynomial f{x, y, z) — f{s, t, u) is 

2 (2a; - 1 + 2s - 2ut^ + 2uH - 2zy^ + 2z'^y) 
{x — s + z^y — zy"^ — u^t + ut^) . 

The first factor fi — 2x ~ 1 + 2s — 2ut^ + 2u^t — 2zy^ + 2z'^y is not symmetric 
near-separated because fi{x,y, z,x,y, z) ^ 0. On the other hand, the second 
factor f2 = x — s + z'^y — zy^ — u^t + ut^ does satisfy f2(x, y, z, x, y, z) — 0. 
Note that by a previous remark, the components of the decomposition can be 
considered as polynomials. Then f^ can be written as fi — h{x, y, z) — h(s, t, u). 
Taking h[x, y, z) ~ f2{x, y, z, 0, 0, 0) = x + z'^y — zy'^, we check that it satisfies the 
previous equation (see theorem\^. The left component g is also a polynomial, 
and by theorem^ has degree 2. Solving the equation f — g{h) we have the 
multi-univariate decomposition: 

(4t2 -2t + 10,x + z'^y - zy"^). 

fn 

Example 3. Let f = ^ with 

fd 

/„ = y2^2 ^ 2x'^yz'^ - 2ySx + z'^x'^ - 2z'^xy^ + 

-81a;2 - AbOxyz - 625y'^z'^, 
fd ^ y^x^ + 2x'^yz'^ - 2y^x + z'^x'^ - 2z'^xy'' + y^" 

-162x2 - 900xyz - 1250y'^z'^. 

We look for all the intermediate fields o/Q(/) C Q{x,y,z) with transcendence 
degree 1 over Q. First, we will factor the polynomial 

fn{x,y,z) fd{s,t,u) - fn{s,t,u) fd{x,y,z) = -625/1/2, 

where 

fi——xtz^u + -^xt^ — zsty — zu^sy + zt^y — ^xz^ s 
— ■^xu'^s — ^xys — xyut — ^xts + ^sy^ + uty^, 

f2=—xtz'^u — -^xt^ + zsty + zu^sy — zt^y — ■^xz'^s 
^-^xu^s — -^xys — xyut + -^xts + -^sy^ + uty^ . 

We have fi{x,y, z,x,y, z) ^ 0, thus fi is not symmetric near-separated. On 
the other hand, f2{x,y,z,x,y,z) — 0. Moreover, 

/2 = -zt^y + uty^ + {^^^^ ^ ~ 2/"* 




9 9 , o 9 , 

— t H y sx. 

25 25 25*^ ' 
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We check that ji is symmetric near-separated, by solving a linear system of 
equations. Define 

f2{x, y, z, 1, 0, 0) = r(x, y, z) = -^xz^ ~ ^xy + ^y^ 

,2b 25 ^ 25 

9 2 9 \ 9 , 

z V X H V 

25 25^y 25^^ 

Next, we compute SQ{y,z) such that 

9 9 

— y5so(t, u) - ^t^so{y, z) = -zf^y + uty^ . 

25 

Let so(j/, z) — 05(2) y^ + ■ ■ ■ + aQ(z). Then ai — — z and oq = 02 = = 

n u 25 J / N ri{y,z)so{t,u) - cio ^ 

04 = 05 U. Hence, sq = — zy ana si{y,z) — = 1. 

9 ro[t,u) 

25 

Thus s — X ^ — 9"'^^' ~ ^ '^^'^ (^1'*) symmetric near-separated 
representation of p: 

9^9 9 c: 

r = xz XV H y 

25,, 25 ^ 25 
25 

s = x + y zy. 

Now we compute g, which is a univariate function with degree 2. Solving the 
corresponding linear system of equations we obtain 

625t^ - 6561 



625t2 - 13122' 



3.3 Second algorithm 

For this algorithm, we suppose that K has sufficiently many elements. If it is 
not the case, then we can decompose / in an algebraic extension M^[lu] of K, and 
then check if it is equivalent to a decomposition with coefficients in IK; this last 
step can be done by solving a system of linear equations in the same fashion 
as the computation of g. The algorithm is based on Corollary [TJ we will need 
several technical results too. 

Lemma 1. Let f G K(x). Then for any admissible monomial ordering > there 
are units u £ Vi G M.{xi), i — l,...,7i such that, if f = fnl fd = 

uo f{vi, . . .,Vn), then Im /„ > Im J^, 7„(0, . . . , 0) = andJ^{Q, ... ,0) ^ 0. 

Proof. Let > be any admissible monomial ordering. Let Ui £ be a unit 

such that /i — fin I /irf = verifies Im /i„ > Im fx^. Such a unit always 

exists: 

- If Im /„ < Im /d, let Ui = 1/y. 
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- If Im /„ - Im /rf, let u, - (1/y) o {y - J^). 

- If Im /„ > Im fd, let ui = y. 

Let a = (ai, . . . ,q;„) G K" such that fid{(£) 7^ (such a a exists if K has 
sufficiently many elements). Let Vi — Xi + ai, i — 1, . . . ,n and /2 = j^nl In = 
. . . , u„). Then /2d(0, . . . , 0) ^ 0, and we can take 

/2»(0,...,0) 

""^ /2d(0,...,0) 

so that f = u o f(vi, . . . ,Vn) verifies all the conditions. □ 

Lemma 2. Let i,j,k £ N with i<j<k. P, QCz K[ac] and > an admissible 
monomial ordering such that Im P > Im Q. Then Im P^Q^^^ > Im P^Q^^^ . 

Lemma 3. Let f = ] ^IJ d ^ K(x) such that Im /„ > Im /,i(0, . . . , 0) = 
and /d(0, ... ,0) 7^ 0. Then, for every uni-multivariate decomposition f — g{h) 
there exists an equivalent decomposition f = g{h) with 'g = 'g^/'g^, deg g„ > 
deg gd and g„(0) = (thus ^^(0) ^ 0). 

Proof. As in the proof of Lemma [H there exists a unit ui such that if hi = 
u{h) = hin/hid, then ft,i„(0, . . . , 0) = 0. Let gi = g{u~^) = (a„?/" + ••• + 
ao)/ibyy'' + ■■■ + bo). Then 



V — U 



and by CorollarylH = (b^h-^^ + ■■■ + boh%)hZr^""'^°^ ._As /^(O, . . . , 0) 
and /ii„(0, . . . , 0) = we must have hid{0, . . . , 0) 7^ 0. But /„(0, . . . , 0) = and 
/„ = {O'uhin + ■ • • + aoh"^)h^^^^^ thus = 0. Next, we will prove that 

there is an equivalent decomposition verifying the condition on the degrees of 
the left-component. To that end, we will consider three cases. Let > be any 
admissible monomial ordering and w — deg gi = max{u,w}. 

— If Im hin < Im hid then using repeatedly Lemma [2l 

Im 7n = Im < Im /li"^ = Im 

which contradicts our hypothesis. 

— If Im hin > Im hid, then applying Lemma [21 

im7„ ^ Im /.^„/.rr^''--''> 

Im/, = lm/iy>7r{"-'''°>. 

As Im /„ > Im fd by hypothesis, by Lemma [2] again we must have u > v, 
that is, deg gin > deg gid. 
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— If Im hin = Im hid then, as in Lemma [TJ we can cancel the leading monomial 
of hid with a unit U2 on the left, so that / — g2{h2) with Im /i2n > Im h2d 
which is the previous case. 

□ 

Let / = g{h) be a uni-multivariate decomposition of / with / = fn/ fd,g = 
{auU^ + • • • + a())/{bvy" + • • • + 60) and h ~ hn/hd- By the previous lemma, we 
can suppose u > v and g{0) — 0, i.e. ao = 0. Then, as 

r ^ auK + ■ ■ ■ + aihnhy^ 
{Khl + ■■■ + bohX~^ 

we have that /i„ | /„ and hd\ fd- This is the key to the following algorithm. 
Algorithm 4. 
Input: / e K(x). 

Output: {g, h) a uni-multivariate decomposition of f , if it exists. 
A Compute u,vi, . . . ,Vn as in Lemma\^ Let 

7 = 7«/7d ^UO f{vi,.. .,Vn) 

B Factor /„ and f ^. Let D = {(Ai, _Bi), . . . , (Am, _Bm)} he the set of pairs 
(A, B) such that A, B divide /„ , fd respectively ( up to product by con- 
stants). Let i — 1. 

C Check if there exists g G IK(y) with f — g{Ai/Bi); if such a g exists, RE- 
TURN {u-\g),h{vi\---,v-^))- 

D If i < m, increase i and go to C, otherwise RETURN NULL (f is uni- 
multivariate indecomposable) . 

Example 4. Let 

f = Az^y'^ - 8z^y^ + Sz'^yx + Az'^y^ - Szy'^x 
+4.t2 - 2z2y + 2zy^ - 2x + 10 

as in Example\M VVe take u = t — 10 G K(t) and vi ^ x,V2 — y, V3 — z. Then 

7 = Az'^y'^ - 8z^y^ + Sz'^yx + Az'^y'^ - Szy'^x 
+4x2 - 2z'^y + 2zy2 _ 2x. 

We factor f = 2{x + z'^y — zy^){2x — 1 + 2z^y — 2zy^). We first take the 
candidate {x + z^y — zy'^). We have to check if there are values of 01,02 for 
which g — 02^^ + ait verifies f = g{x + z'^y — zy"^). We find the solution 
a2 = 4, 01 = -2. Thus f ^ {At^ - 2t + I0){x + z^y - zy^). 



14 



Example 5. Let f ~ — with 

Id 

fn = y^x^ + 2x'^yz'^ - 2y^x + z'^x^ - 2z'^xy^ 

- 81a;2 - 450xyz - &2by'^z'^, 

fd = y^x^ + 2x^yz^ - 2y^x + z'^x^ - 2z'^xy^ 

- 162a;2 - 900xyz - 1250y^z^, 

as in Example \3i Let > be the pure lexicographical ordering with y > x > z. 
Then Im /„ = Im = j/^". Following the proof of lemma[]\ let ui — l/(t — 1), 
then ui{f) = fin/ fid with 

/i„ = y2^2 _^ 2x^yz^ - 2y^x + z'^x^ - 2z^xy^ 

+2;^° - 162x2 - gOOxyz - 12502;2z^ 
/irf = 81a;2 +450a;?/z + 625?/2z2. 

Now, let a = (1,0,0), so that the denominator of the previous expression is 
non-zero at the point a. Then f2n/f2d — u{ f{x + 1, j/, z)) with 

f2n — y^x'^ + 2y'^x + y2 + 2x'^yz'^ + -iyz'^x + 2yz'^ 
-2y^x - 2y^ + z'^x^ + 2z'^x + z^ ~ 2z'^xy^ 
-2z^y^ + - 162x2 - 324x - 162 
-mOxyz - 900yz - UbOy'^z^, 

f2d = 81x2 + 162x + 81 + 450xj/z + 450yz + 625j/2z2. 

As /2„(0, 0, 0) = -162 and /2d(0, 0, 0) = 81, ifu2 = t + 2, we have that 
U2{uiif{x + l,y,z)))=7=^ 

fd 

verifies the conditions of Lemma\^ We factor /„ and f ^: 

7_n = {z^ + z'^x + y + xy ~ y^f , 
fa ^{9x + 9 + 25yzf . 

As the degree is multiplicative and deg / = 10, and also Im /i„ > Im hd, the 
possible values of hmhd are 

hn = z^ + z^x + y + xy — y^, 

hd e {l,9x + 9 + 25yz,(9x + 9 + 25yz)^} 

To check them, let g = . We substitute h ing and solve the homogeneous 

linear system obtained by comparing the coefficients with those of f . 

" If hd — 1, there is only the trivial solution, thus h is not a candidate for f. 

— Lf hd = 9x + 9 + 25yz, we get the non-trivial solution 
a2 = bo — I, fli = &i = 0, thus f has a uni-multivariate decomposition 

iu^\u^\g)) , h{x - l,y,z)) = (^, 'tiff/ ). 



15 



Table 1: Average computing times (in seconds) 



n 


d 


AlgEl 


Fact. 


AlgH 


Fact. 


2 


10 


32.17 


23.15 


27.03 


22.44 


2 


25 


68.20 


46.34 


51.10 


40.33 


2 


30 


89.40 


62.48 


91.22 


71.06 


4 


8 


54.37 


38.56 


32.07 


25.47 


4 


25 


89.75 


65.95 


64.41 


46.72 


4 


30 


156.87 


110.30 


134.60 


99.87 


8 


10 


234.90 


162.89 


156.12 


116.66 


8 


25 


349.44 


235.41 


341.11 


276.85 


8 


30 


654.72 


454.36 


678.89 


511.01 



^ If hd = (9a; + 9 + 25yz) , the only solution is the trivial one. 

Therefore, any uni-multivariate decomposition of f is equivalent to the decom- 
position [g, h) computed before. 

4 Performance 

Both algorithms run in exponential time, since the number of candidates to be 
tested is, in the worst case, exponential in the degree of the input; the rest of the 
steps in both algorithms are in polynomial time. However, in practical examples 
it seems that most of the time is spent in the factorization of the associated 
symmetric near-separated polynomial, in Algorithm [3l or the numerator and 
denominator in Algorithm [H We show in Table 1 the average times obtained 
by running implementations of these algorithms in Maple VI (see [3]) on a 
collection of random functions. The parameters are the number of variables n 
and the degree of the rational function d. We have also included the factorization 
times for each algorithm. 
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